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Abstract: We present an extension of the antenna subtraction formalism at NLO to 
include massive final state fermions. The basic ingredients to the subtraction terms, the 
NLO massive final-final antenna functions are derived and integrated over the correspond- 
ing factorised phase space. Those antenna functions account for all soft, collinear and 
quasi-collinear limits of the QCD matrix elements involving massive fermions in the final 
state. 
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1. Introduction 

Hard scattering processes leading to final states with heavy particles are important observ- 
ables for present and future colliders. Reliable theoretical predictions for these observables 
require the calculation of at least the next-to- leading order QCD corrections. 

For hard scattering processes, involving massless or massive QCD partons, the pertur- 
bative corrections to a given process and at a given order in QCD are obtained when all 
partonic channels contributing to that order are summed. In general, each partonic chan- 
nel contains both ultraviolet and infrared (soft and collinear) divergences. The ultraviolet 
poles are removed by renormalisation in each channel. The remaining soft and collinear 
poles cancel among each other when all partonic contributions are summed over [1]. At 
NLO, one has to combine virtual-one loop calculations with real emission contributions 
from unresolved partons. While infrared singularities from purely virtual corrections are 
obtained immediately after integration over the loop momenta, their extraction is more 
involved for the real emission contributions. Here, the infrared singularities only become 
explicit after integrating the real radiation matrix elements over the phase space. To build 
a general-purpose Monte Carlo program to evaluate observables at NLO requires therefore 
an analytic cancellation of infrared singularities before any numerical integration can be 
performed. 

For the task of NLO calculations, several systematic and process-independent proce- 
dures are available and have been applied to a variety of processes [2]. The two main 
methods are phase space slicing [3, 4] and the subtraction methods [5, 6]. Except for the 
phase space slicing technique, all subtraction methods consist in introducing terms which 
are subtracted from the real radiation part at each phase space point. These subtrac- 
tion terms approximate the matrix element in all singular limits and should be sufficiently 
simple to be integrated over the corresponding phase space part analytically. After this 
integration, the infrared divergences of the subtraction terms become explicit and the inte- 
grated subtraction terms can be added to the virtual corrections, thus yielding an infrared 
finite result. 

One of the widely used subtraction formalisms at NLO is the dipole subtraction for- 
malism of Catani and Seymour, which in its original formulation [6] deals with massless 
partons in final and/or initial state at NLO. 

Another subtraction scheme is the antenna formalism [7, 8, 9], which constructs the 
subtraction terms from so-called antenna functions. The antenna functions describe all 
unresolved partonic (soft and collinear) radiation between a hard pair of colour-ordered 
partons, the radiators. These functions can be derived systematically from physical matrix 
elements, as shown in [10, 11] and can be integrated over the factorised phase space. In 
the antenna subtraction method, the subtraction terms are constructed from products of 
antenna functions with reduced matrix elements (with fewer partons than the original 
matrix element) such that these subtraction terms are acting on individual colour-ordered 
real radiation matrix elements. So far, this formalism can handle massless partons in final 
and/or initial state at NLO [7, 9] and it can treat massless final state parton radiation at 
NNLO [8]. It is the only formalism where the NNLO formulation has been worked out in 
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full and applied to evaluate NNLO corrections to jet observables, namely for e + e — ► 3 
jets [12, 13]. 

To evaluate hard scattering processes involving massive hard partons in the final state 
at the next-to-leading order level, a subtraction formalism is needed as well in order to 
isolate and cancel infrared divergences among different partonic contributions to the cross 
section. However the extension from massless partons to massive partons involve compli- 
cations and new features. The kinematics is more involved due to the finite value of the 
parton masses. 

Furthermore, QCD radiation from massive partons, although not leading to strict 
collinear divergences, (which are regulated by the mass of the massive particle), will be 
proportional to In Q 2 /M 2 , where M is the parton mass and Q is the typical scale of the 
hard scattering process. In kinematical configurations where Q 3> M, these logarithmi- 
cally enhanced contributions become large and can spoil the numerical convergence of the 
calculation. Although these logarithmic terms cancel in the final NLO result they appear 
at intermediate steps of the calculation. Those can be easily traced though. Indeed, these 
logarithmic enhanced terms are related to a process-independent behaviour of the matrix 
elements; its singular behaviour in the massless limit (M — > 0). This singular behaviour is 
related to the quasi- collinear [14] limit of the matrix element. For massless partons, infrared 
divergences arise in soft and collinear kinematical configurations and in these configura- 
tions matrix element and phase space obey QCD factorisation formulae [15]. Similarly, for 
massive partons, matrix element and phase space will obey factorisation formulae in the 
soft and quasi-collinear limit. 

The dipole subtraction formalism of Catani and Seymour has been extended to include 
effects of massive partons in [14, 16, 17] at NLO. An extension of the antenna subtraction 
method including the radiation effects of massive partons has been so far missing. 

It is the purpose of this paper to present an extension of the antenna subtraction 
method to include radiation of final state massive fermions in order to be able to evaluate 
the production of massive particles beyond leading order within this formalism. For the 
sake of clarity, in this paper we restrict ourselves to the kinematical situation of a colour- 
neutral particle decaying into one or two coloured massive fermions of equal masses. In this 
situation, infrared singularities appear only due to final state radiation and only two scales 
are involved in the problem. As the radiators are in the final state only, the corresponding 
antenna functions involving those are regarded as final-final antennae. 

The basic structure of the extension of the antenna subtraction formalism to include 
massive particles can be carried over to partons in the initial state. However, to deal 
with the production of massive final states at hadron colliders, the construction of the 
appropriate subtraction terms and their integrations cannot be solely performed with the 
ingredients presented in this paper. Besides massless initial-initial antenna functions de- 
rived in [9], which involve massless initial state radiators, and massive final-final antenna 
functions which will be derived in this paper, antenna functions involving a massless initial 
state and a massive final state are also required. As in the massless case, those initial- 
final antenna functions can be obtained from the final-final massive antenna functions by 
crossing one massless parton from the final to the initial state. Their integration over the 
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corresponding phase space cannot be taken over from the final-final case, though. 

Furthermore, the extension of the antenna subtraction method involving massive final 
state particles is not restricted to Standard Model processes. It can also be applied to the 
production of supersymmetric particles, for example. An appropriate treatment of massive 
final state scalars, like squarks could also be realised within this framework. 

So far, final-final massless antennae have found some interesting applications: Starting 
from the three-parton antenna functions given in [8], a parton shower event generator 
VINCIA [18] has been constructed. One could therefore envisage a further implementation 
in VINCIA of the massive NLO final-final antennae presented in this paper. A combination 
between NLO calculation in the antenna subtraction formalism and parton showers could 
then become feasible for processes involving massless and massive final state fermions. 
Originally, parton showers based on the antennae approach for final state particles had been 
derived for the event generator ARIADNE [19]. Furthermore, parton showers related to the 
dipole subtraction formalism involving dipole functions were considered in the litterature 
[20, 21, 22] and are already implemented [23] in SHERPA [24]. 

The paper will be organised as follows. In Section 2, we recall the structure of the 
antenna subtraction terms at NLO for massless final state partons and state how it extends 
with the presence of massive partons. Section 3 presents the formulae for the massive 
antenna functions while in Section 4 a list of all non-vanishing limits for these antennae 
is given. In Section 5 we present the results for the integrated massive antenna functions. 
To validate our results for the integrated massive quark-antiquark antenna function, we 
recompute the NLO correction to the hadronic decay rate in Section 6. Section 7 contains 
our conclusions. 

2. Antenna subtraction with massive final-final configurations 

In configurations involving final-final antennae, both radiators are in the final state. For 
massless radiators, this case was described in detail at NLO in [7] and NNLO in [8]. For 
massive radiators, the same factorisation formulae for the subtraction terms will hold. Only 
the ingredients, phase space and antenna functions will be modified to take into account 
the mass of radiators. For clarity reasons, we will keep the same notation as in [8]. 

The subtraction term for an unresolved parton j, emitted between massless or massive 
hard final-state radiators i and k reads, 



This subtraction term involves the m parton amplitude A4 m which depends only on the 
redefined on-shell momenta k±, ki, kx, ■■, ^m+i- The momenta kj,kx are linear combi- 
nations of the momenta ki , kj , kk while the tree- level antenna function Xf- k depends only 




Xijk \Mmih, • • • , Kj, K K , . . . , k m+1 ; q)\ 2 J^\k u . . . , K u K K , . . . , k m+1 ;q) . 



j 



(2.1) 
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on ki,kj,kk, X stands for different antenna types denoted (in the massive NLO case) by 
A, E, D or G depending on which radiators are involved. The jet function Jm in eq.(2.1) 
does not depend on the individual momenta hi, kj and k^, but only on Ki,Kk- This 
function ensures that any partonic contribution give rise to an m-jet final state. 
The phase space d<l> m+ i can be factorised as follows, 

d$ m+ i(fci, . . .,k m+ i;q) = 
d$ m (k 1 ,...,K I ,K K ,...,k m +i;q) ■ d^x^iki^^k^Kj + K K ) . (2.2) 

For massive radiators, the phase space d$ m is the ci-dimensional {d = 4 — 2e) phase space 
for m outgoing particles with momenta pi, - ■ • ,p m , masses mi, • • • ,m m with total four- 
momentum q^. It reads, 

d *™(*>-->^*>=2iSN ■■■ dS^^ *'<«-*'-■■■-*">• (2 - 3> 

with 



Ei = J\pi\ 2 + ml (2.4) 



S m is a symmetry factor for identical partons in the final state. d&x ijk is the NLO antenna 
phase space, it is proportional to the three-particle phase space relevant to a 1 — ► 3 decay. 
This can be seen by using m = 2 in the above formula (2.2) and exploiting the fact that 
the two-particle phase space is a constant. The massive NLO antenna phase space will be 
derived in Section 5. 

The tree-level antenna function X^- k describes all the situations where parton j is 
unresolved. Those are obtained by normalising the corresponding colour-ordered three- 
parton tree-level squared matrix elements to the squared matrix element for the basic 
two-parton process, omitting all couplings and colour factors. Indeed, the normalisation 
for the antennae is chosen such that in the unresolved limits those yield exactly the well- 
known collinear splitting functions and soft eikonal factors. As such the massless and 
massive antenna functions have mass dimension —2 and are defined by [8], 

|_A/fO fc |2 

X ijk = SijhJK i u0 | 2 , ( 2 -5) 
\ JVI IK\ 

where S denotes the symmetry factor associated to the antenna, which accounts both for 
potential identical particle symmetries and for the presence of more than one antenna in 
the basic two-parton process. The massive antenna functions will be given below in Section 
3 while their unresolved limits will be presented in Section 4. 



3. Massive NLO antenna functions 

The massive final-final antenna functions can be derived from physical matrix elements 
of the same processes as in the massless case [8, 10, 11] but considering the final state 
radiators to be massive. More precisely, the quark-antiquark antenna function which we 
denote by A° Q q (to be distinguished from A® g) can be derived from 7* — ► QQ + g with Q 
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being a massive quark of mass itiq, g is a massless gluon. We find three types of massive 
quark-gluon antennae involving either one or two massive partons. Those can be derived 
from the following processes: (x — ► g+ 2 partons). For these final state partons being two 
gluons, we derive Dg gg while for these partons being a quark-antiquark pair, we derive 
EQ q ,q, and E®q,q,. Lastly, we can also define a massive gluon-gluon antenna through the 
process H — > gQQ, with one of the gluons in the final state emitting a massive QQ pair. 
All new massive antenna functions are given below with mg denoting the mass of the 
massive quark Q, g is a massless gluon and the invariant is chosen to be Sij = 2 p, L • pj 
making the mass dependence in the formulae below explicit. E cm is always the rest energy 
of the decaying particle. All antennae listed below are known exactly at all orders in e. 
The O (e) parts have been omitted for conciseness only, those are however already needed 
at this next-to-leading order level. 

The quark-antiquark massive antenna function Aq^q reads, 



^, (1,3,2) = _ T 4(^L. + — + — + - + 

V an J 

2 / 8S12 2si2 2si2 2s 2 3 2 2 2s i3 \ C3 

+ m Q I 2 2 2 T~ 

\ S 13S23 «i 3 «23 s 13 s 13 s 23 «23 



2s? 2 2S12 , 2si2 , S 2 3 , «13 



+ ^(-i-i)) +0(£) ' 



This function has been normalised to the two-particle matrix element relevant for 7* — > QQ, 
whose matrix element squared (omitting couplings) is given by 



A °QQ^ 2 ) = 4 [(! " e ) ^ + 2 m^] . (3.2) 



The quark-gluon massive antennae with either two gluons or a massless quark-antiquark 
pair in the final state are normalised by the two-particle matrix element squared relevant 
for x 99-, w ith the gluino g being massive with mass tuq, the gluon g being massless. 
This two-particle matrix element squared omitting couplings reads: 



*yi,2)=4(l- € ) (E 2 cm -m 2 Q ) 2 , (3.3) 



where X can stand for E or D here. The quark-gluon massive antenna Dg gg with two 
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gluons and a massive radiator Q in the final state reads, 



n , 1 ( ( 4s?n Ash 4s?, 2s?, 

D° Qgg 1,3,4 =77^2 2T2 4 9 Sl3 + 9 Sl4 + ^ + ^ + ^ + ^^- 



6S13S14 , 4sf 4 2sf 4 3si 3 s 34 3s 14 s 34 s§ 4 s§ 4 
H 1 h os 34 H 1 1 h 



S34 «34 Si 3 S 34 S14 S13 S13 S14 

• i , 2s i3 , 4s i3 , 4 «i4 , 2sf 4 6s 34 4si 3 s 34 
-m ( 6 + — 2^ H 1 h -2— H 1 5 

Sf 4 S14 S13 Sf 3 S13 Sf 4 



6S34 4si 4 s 34 2s|4 2s|4 2s| 4 \ 

S14 Sf 3 Sf 3 Sf 4 Si 3 S 14: J 

+2E cm m Q - E cm m% 2534 + mi—) + O (e) . 

v Sl 3 Sl 4 V S1 3 S14/ 



(3.4) 



The quark-gluon massive antenna Eq q ,-, with a pair of massless quarks and a massive 
radiator quark Q in the final state reads, 

1) - ^L^ 4 (l. + •« + J + £ - + O («) . (3.5) 

In the case where the gluino is massless and the gluon radiates a massive quark- 
antiquark pair, the two-particle matrix element which serves to normalise the quark-gluon 
massive antenna E°„,=, is different. It reads: 

qQ'Q' 

E qg (l,2)=4Ei n (l-e). (3.6) 

The quark-gluon massive antenna E®q,q, with a pair of massive quarks and a massless 
quark-radiator in the final state reads: 

Kq'Q'^ 3 ' 4 ) = 4iL (, 34 + 2mlr (S?3S34 + S?4 " 34 + Sl3 " 34 + Sl4 ' 34 

(3 7) 

+m 2 Q (4s^3 + 4si 3 si 4 + 4s^ 4 + 6si 3 s 34 + 6si 4 s 34 ) 
+8m 4 Q (ai 3 + si 4 )) + O (e) . 

Finally, the gluon-gluon massive antenna where the final state gluon emits a massive quark- 
antiquark pair reads: 

GgQ'Q'C 1 ' 3 ' 4 ) = 4^4 ( S34 + 4 2m 2)2 ( S ?3 S 34 + S? 4 «34 + 4m^ (s? 3 + s\ A + S13S14)) + O (e) . 

Q (3.8) 
This antenna function has been normalised to the two-particle matrix element relevant for 
H ^ 99 whose matrix element squared is given by, 

G 9 ,(l,2)=4£ c 4 m (l-e). (3.9) 
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4. Singular limits 



The antenna functions listed above encapsulate all single unresolved limits of tree-level 
QCD matrix elements involving massive final states. The factorisation properties of tree- 
level QCD squared matrix elements are well-known [15] for massless and for massive final 
state partons. In the antenna picture, for massless radiators, the unresolved and massless 
parton emitted between those can be either collinear or soft. In these corresponding unre- 
solved limits, defined as in [8] the (m + l)-parton matrix element squared factorises into 
a reduced m-parton matrix element and the well known soft eikonal factor when a gluon 
is soft and one of the three different Altarelli-Parisi splitting functions [25] when two of 
the three massless partons are collinear. When the radiators which emit a massless and 
unresolved parton between them are massive instead, the limiting behaviour of the matrix 
elements has to be changed and the unresolved factors have to be modified accordingly. In 
this section we first recall the single unresolved massless factors and present their exten- 
sion needed to deal with massive radiators. Finally, we list all non-vanishing limits of the 
three-parton massive antennae defined in Section 3. 

4.1 Single unresolved massless factors 

When a soft gluon (j) is emitted between two hard and massless partons (i and k), the 
eikonal factor Sijk factorises off the squared matrix element, 

Sak = ^- (4-1) 

When two massless partons become collinear, we have different splitting functions corre- 
sponding to various final state configurations: a massless quark splits into a quark and a 
gluon (P ?s _>q), a gluon splits into a quark-antiquark pair (P q q^a) or a gluon splits into 
two gluons {Pgg^c)- These are given by: 



p qg~*Q( z ) 



1 + (1 - z) 2 - ez 2 

z 

z 2 + (l-z) 2 -e 



Pgg^a(z) = 2 + L-i + z(l ~ Z)j . (4.2) 

In these equations, z is the momentum fraction carried by the unresolved parton and the 
label q present in these splitting functions can stand for a quark or an antiquark: P qg ^Q = 
Pqg^Q by charge conjugation. Q or G appearing in the collinear splitting functions denotes 
the parent particle of the two collinear partons % and j. 

4.2 Single unresolved massive factors 

When massive partons are present in the final state, those can be soft but they cannot be 
strictly collinear, the mass is regularizing the collinear divergence. The relation between 
the massive matrix element squared and the splitting functions needs to be extended from 
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massless to massive. Similar factorisation formulae as in the massless case will hold pro- 
vided the collinear limit is generalized [14] to the quasi- collinear limit. 

The limit when a massive parton P of momentum p^ decays quasi-collinearly into two 
massive partons j and k is defined by, 



p^zp^p^(l-z)p^ (4.3) 
"*(?*)■ 



P 2 = ml k) . (4.4) 



with the constraints, 
at fixed ratios, 



Pj • Pk,mj,m k ,m jk -> (4.5) 



j3_ jsL (46) 

Pj ■ Pk Pj ■ Pk Pj ■ Pk 

The key difference between the massless collinear limit and the quasi-collinear limit is 
given by the constraint that the on-shell masses squared have to be kept of the same order 
as the invariant mass (pj +Pk) 2 , with the latter becoming small. In this corresponding 
quasi-collinear limits, the (m + l)-parton matrix element squared factorises into a reduced 
m-parton matrix element and unresolved massive factors. 

More precisely, the single unresolved massive factors in which the real matrix elements 
factorise in the soft and quasi-collinear limits are generalizations of the massless soft and 
collinear Altarelli-Parisi splitting functions defined above. The massive splitting functions 
denoted by Pij^^(z, /J, 2 j) for parton (ij) splitting into partons i and j in d-dimensions 
(d = 4 — 2e) will be given below. In these expressions, z is the momentum fraction carried 
by the unresolved parton becoming quasi-collinear to the massive parton. All the mass 
dependence can be parameterized by fijk = (m 2 + rn\)/[{pj + Pk) 2 — m (jk)\- ^ ne number 
of gluon polarizations is chosen within the dimensional regularisation procedure as d — 2 . 
Those generalized massive splitting functions were also given in the appendix of [16] and 
read, 



.2 . / 1 +(■' - -) L - ' - ,,,.2 ) 

m 



P qg ^ Q (z,H qg ) = 2(1q 



Pqq-+G(z, figg) 



2 

2 2 " l Q 

where fi qg is defined by n qg = 

( z 2 + (l-z) 2 -e- 

2 2 2m Q 

where \i qq is defined by [i qq - = 



s q q + 2m 2 Q ] 

The gluon-gluon splitting function P gg ^G( z ) 1S left unchanged. The generalized soft eikonal 
factor Sijk(rrii, m^) depends on the invariants si m = 2pi -p m built with the partons i, j and 
k but also on the masses of partons i and k, mi and nik respectively. It is given by, 

c ( \ 2sik 2m i 2m k (A 7\ 

S ijk (mi,m k ) = ^- (4.7) 

SijSjk s- k 



-8- 



4.3 Singular limits of the massive NLO antenna functions 

We list here the non- vanishing soft, collinear and quasi-collinear limits of the massive 
antenna functions given in Section 3. The singular limits of A^^axe 

A QgQ (!> 3 ' 2 ) — ^ 5 132("T-Q, m Q ), 

A Q 9 Q (1. 3 ' 2 ) ^ ^ (48) 

^Q fl Q (!. 3 ' 2 ) ^ ^ 

For ^Qss^ have 

^(1,3,4)^^34^,0), 
£>S M (1,3,4) ^ ^-P,^ (z,^ ff ), 

2?0 M (l,3,4)i^— 

The only non- vanishing singular limit of Pg^^is the collinear massless limit of the 
massless quark-antiquark pair, 

E° Qqlq , (1,3,4) -Lp Q ^ G ( z ), (4.10) 

s 34 



(4.9) 



while the only non- vanishing limit of E®q,q, is, 



Kq>Q> ^ 3 ' 4 ) (4-11) 
Finally G^^does also have a quasi-collinear quark-antiquark limit given by 

G° gQ/Q/ (1,3,4) S3 J 2m2 P^G(z,^). (4.12) 

Q 

5. Integrated massive NLO antenna functions 

In this section we present the results for the integration over the appropriate antenna phase 
space of the massive final- final three-parton antenna functions defined in Section 3. 

As in the massless case, the phase space associated to the subtraction term given in 
eq.(2.2) can be factorised. In a first step, we present the phase space factorisation of a 
three-particle phase space with three massive particles into a two-particle phase space and 
an antenna phase space. The latter is proportional to a three particle massive phase space. 
We need to consider two cases: three particle final state with one massless parton and 
either two massive partons of equal mass, or one massive parton and two massless partons. 
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5.1 Phase space factorisation 

In the rest frame of a decaying photon of momentum q, the <i-dimensional two-particle phase 
space d&2(pi,PK', q) given in eq.(2.3) for two outgoing particles / and K with momenta pi 
and pk and masses mi and mx reads: 

d* 2 ((pi,PK; q) = (2n) 2 - d (£) 0O 1_d/2 ((E 2 cm - m 2 - m 2 K f - Amjm 2 ,) ^ <H2 d _i, 

(5.1) 

with 

2ir d / 2 

dfi -w (5 - 2) 

where we have used 



/ 



/ ^2ET = \\ dE i\Pi\ d - 3 <M d -i. (5.3) 

In this last equation, f2<j_i parameterizes the solid angle of the d — 1 components of pi in 
d — 1 dimensions. The three particle phase space d&3(pi,pj,p k ;q) may be written as the 
product of a two particle phase space d$2(pi,PK', q) defined above and an antenna phase 
space d^x ijk which depends on the momenta Pi,Pj and p k and on the masses mi,mj and 
m k . 

Using spherical coordinates in d dimensions, we have 

f = I ^[b1||p>in<- 3 d^di^dfldJWJVi (5.4) 

with being the angle between the vectors pi and pj. fid-i> parameterize the solid 

angles of the d — 1 and <i — 2 components of the vectors and pj respectively, in d — 1 
dimensions. We find, 

J d<S> z (jpi,pj,p k -q) = j d<f> 2 (pi,PK;q) x (5.5) 
with the antenna phase space d&x ijk given by, 



/ 



d$x ijk (sij,s jk ,s ik ) = 

{2,r d ^-\ ((E 2 cm - m] - m 2 K ) 2 - Am 2 m 2 K ) " 

y dsij ds jk ds ik S(E 2 m - m 2 -m 2 -m\- s„ - s jfe - s ik ) 

[4A 3 (p i ,p j ,p k )]— 6(A 3 (p i ,p j ,p k )). 

The masses mi and appearing in this equation are combinations of the masses mi, mj 
and m k and the integral over Q.d-2 has been performed. The function As(pi,Pj,p k ) is the 
Gram determinant for massive particles of momenta Pi,Pj,p k given in terms of invariants 
Sij = 2pi ■ pj and masses m,, mj,m k by, 

^(Pi,Pj,Pk) = t (sijS ik s jk - m 2 s 2 k - m\s 2 j - m 2 s 2 k + 4m 2 m 2 m 2 k ) (5.7) 
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As in the massless case, the factorisation of the phase space is exact. The expression 
for the antenna phase space d$x ijk given in eq. (5.6) can also be found in [14, 17], where 
it was called a dipole phase space. At NLO, both massless and massive dipole and antenna 
phase spaces are the same. The NLO massive antenna functions differ however from the 
NLO massive dipoles functions given [14] as did the NLO massless antenna [8] from the 
massless dipole functions [6]. 

In the following, we will only consider two kinematical configurations of the final 
state particles. For the case where one of mass vanishes (m,j = 0) and the two other 
masses are equal (rrij = mt = itiq), i.e mj = uik = rrii in eq.(5.6), we use the following 
parameterization, 



J drr 1_2e (l - r) _e+ 3 (1 - r r)~^ J ds (l - s 2 )' 



1 4 ™Q 
r = l--^-, 

cm 

s ij = 2 E c™ r ° r 



s 3k = 2 E cm r r 




(5.8) 



and for the case where two particles are massless m; = rrij = and the third one is massive, 
m k = "m-Qi m which case mi = mj = 0, tuk = in eq.(5.6), we use, 



.*(0,0,m) _ {^Y (F 2 f.. n )2-2e 



f 1 duu 1 - 2 ^! - uy- 2e {l - u u)- 1+e f 1 dvv- e (l - 
Jo Jo 



vy 



m 2 Q (5.9) 

cm 

2 2 U(1-U)V 
V ~ Cm ° (1 - UQU) ' 

Sik = E 2 m u {l - u). 

5.2 Integrals 

Following the extension of the integration by parts method [26, 27] in [28, 29], to reduce 
the number of real phase space integrals, we have expressed all invariants in terms of 
massive propagators and expressed the three on-shell conditions pf = mf, i = 1,2,3 as 
cut propagators. The reduction to master integrals using the Laporta algorithm [30] was 
done independently, once with an in-house implementation in FORM [31] and once with 
the mathematica package FIRE [32]. We find five master integrals, four of which can be 
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evaluated in terms of hypergeometric functions for arbitrary e. For the last one, only the 
expansion up to order e 2 will be given. Differential equation techniques [33, 34] are used 
to compute this master integral. 

The integrated antennae can be separated in two categories corresponding to the phase 
space parameterizations used. The integrated antennae with two massive particles are 
obtained by integrating the antennae with two massive final state partons presented in 
Section 3 over the antenna phase space given in eq.(5.8). Those can be expressed in terms 
of two master integrals /j m '°' m ) (i=l,2) and are given exactly at any order in e by, 



A 



£ c 2 m e(l-2e)r (l-r )(3-r -2e) 
((2(1 - e)(-15 + 8r - r 2 + e(28 - 12r + 4r 2 ) - e 2 (12 + 4r + 4r 2 ) + 8e 3 r )) l[ mfi ' m) 



+ -|-24(1 - e) 2 (3 - r - 6(2 + r ) + 2 e 2 r )/f '°' m) 



£° - - — 
qQ'Q' ~ E 2 



1 /2ro(3 + ro — 4ero) r ( m ,o,m) 1 ,(l+?"o — 2ero) j-(m,0,m) \ 
^ { W^o) 1 ~ W m (1 - ro) 2 J ' 

C ft _±_ ( 8 U ~ e ) r l T (m,0,m) 1 n _ ^ rp T (m,0,m)\ 

W - El m \ 3(1 - r ) ^ E 2 m 8(1 e) (l-r ) 72 )' 



(5.10) 



The two master integrals /^ m '°' m ) are given exactly below. _r( m >°> m ) j s the integrated phase 
space measure and is given by 

j(rn,0,m) _ / ^(m,0,m) 

_ (F 2 q-.2-2.g-2. -2 + . r(2-2e)r(3-3e) (1 7 \ 

(5.11) 



and 



_ rF 2 , 2 -e r 3-2 e oi-2. -2 +e r(3 - 2e)r(4 - 3e) A 9 
-(*U ^o 2 vr r(8 _ 6e) 2 Fi^-,3-26,--3e,r 



(5.12) 



with p 2 = m 2 and p 2 = 0. 

The integrated antennae with one massive particle are obtained by integrating the 
antenna functions with one massive final state presented in Section 3 with the phase space 
given m eq. (5.9) (with ji = Those can be expressed in term of three master integrals 
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jO,0,m ( i=lr-3 ) ag follows, 

V ° Q " = E^THl-e)(l-2eW(l-^(l + ^ + ^ 

+2e(l + n)(l - 6/ - 28/) - e 2 (7 + 7/x - 42/ - 40/ - 111/ - 125/) 
+2e 3 (5 + 3/x - 31/ - 22/ - 44/ - 59/) 
-e 4 (l + /i)(7 - 8/i - 38/ + 30/ - 19/) 
+2e 5 (l + / 2 (1 - 3/x - / + 9a* 3 )) /f'°' m) 

+ EL^ (l-2 g )/(l-/ 2 (l + /3 (12/d + / - 2,(1 + M )(3 + 22/) (5-13) 
+e 2 (15 + 15/i + 65/ + 77/) - e 3 (l + /(15 - 12 ft + 3l/) 
+26 4 (l + /(3-6^-7/))4 0Am) 

- ^t^ttt; (/! + 2 / + ^ ~ ^ " ^ + / 2 ) 4°'°' m) , 

cO 2 ( ( <2u 2 + fClu 2 uW 7" ( °'°' m) 2 ~ 3e /WmA 

The master integral jj ' '" 1 ) corresponding to the integrated phase space and the master 
integral J^ ' '" 1 ) are given exactly by, 

/ (0,0,m) = f d$ (0,0,m) 

= (^L) 1 - e ^- 2e 2- 4+2e 7r- 2+e r(2 6) 2* (1 - e, 2 - &, 4 - 4 £ ; .o) , 

(5.14) 



r (0,0,m) 
J 2 



= {El m ?-^2-^-^ T{2 r ^ e) 6) 2 F, (1 - e, 2 - 2e, 5 - 4e; u ) , 

(5.15) 

with pf = 0,p1 = m\. 

The third master integral I^'°' m ^ is defined as follows, 



7 (o,o,m) = f 1 ^(0,0 /- 



kSjk 

(S 2 m )- 1 - e (47r)- 2+e Uo 2e / / , , u- 2e (l-n)- 2e (l-u «)^- e (l-«)- 



dudv 



T(l - e) Jo Jo 1 - ""o^ - (1 - u)u t> 

(5.16) 



To get the integrated quark-gluon antenna denoted by £>g gg up to finite terms, only the 
constant term of the integral /(°' ' m ) i s needed. The evaluation of this integral yields hyper- 
geometric functions of argument uq. The expansion around e = of these hypergeometric 
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functions leads to one-dimensional harmonic polylogarithms [35, 36, 37], H PL-functions, 
denoted here by H. The constant term in /^°'°' m ) yields, 



r (0,0,m) 
J 3 



*-^iftk^< 1 - 1 :"»» + ff ft"»»- (5 ' 17) 



Foreseeing that the results we obtained for the integrated NLO massive antenna func- 
tions could be used as an input for a further development of NNLO antenna subtraction 
involving massive radiators we derived the order e and e 2 of this master integral /(°' ' m ). 

For this purpose, we derived differential equations for I3 with respect to E 2 m and uq. 
The equations obtained are homogeneous with respect to E 2 m and inhomogeneous with 
respect to uq. In the differential equation for ^3 with respect to no the master integrals l\ 
and I2 appear as coefficients of the inhomogeneous part of the equation. By inserting their 
known results, we could solve the differential equation for 23 order by order in e. Some 
care has to be taken however. These master integrals and I3 are obtained after an 

integration over the antenna phase space d^^'°'^ . The differential equation with respect 
to uq however concerns integrals over the full three-parton phase space. 

As we saw in eq. (5.5) the latter phase space can be written as the product of the 
antenna phase space and a two-particle phase space. However, since the two particle phase 
space depends on uq, special care has to be taken solving the differential equation for uq. 
Solving this differential equation will yield the wanted master integral ^3 multiplied by the 
two-particle phase space with one massive and one massless particle denoted by ^(m, 0). 
The master integrals I\ and I2 defined above will not appear as such in the differential 
equation but those appear as multiplied by the two-particle phase space $2(^0). 

More precisely, the two-particle phase space with one massive and one massless particle 
reads, 

$ 2 (m, 0) = 2 -3 + 2e ( ^ r l +£ r(l-^ K)1 „ 2£ ( ^ 2j - e _ (51g) 

The master integrals I\ (uo) and ^2(^0) appearing in the differential equation for I3 are 
obtained as the product of the original master integrals I\ and I2 with the two particle 
phase space measure $2(^,0). Those take the form, 

/i(«o) = O cc {El m ) l ~ 2e (uo) 3-4e 2 i ? i(l-e,2-2e,4-4e;n ) 
h(u ) = O cc (E 2 cm f- 2e (n ) 4 " 4e 2 i ? i(l-e,2-2e,5-4e;n ) 

(5.19) 

with the overall normalisation factor O cc given by, 

O cc = 2^(vr)^^i|. (5.20) 
With this notation, the master integral ^3 is given up to order e 2 by, 

'' = W) oMr " 2 ' (5 - 21> 
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with, 



I = (6H (2; n ) + 6H (1, 1; n ) - 2e (16H (2; u ) + 6H (3; n ) + 16H (1, 1; n ) 

-3# (1, 2; n ) + 12H (2, 0; n ) + 3H (2, 1; u ) + 12H (1, 1, 0; n ) - 6# (1, 1, 1; n )) 

+2e 2 (16ff (2; n ) + 32H (3; n ) + 12H (4; n ) + 16tf (1, 1; u ) - 16// (1, 2; n ) 

-6ff (1, 3; n ) + UH (2, 0; n ) + 16# (2, 1; n ) + 24H (3, 0; n ) + 6H (3, 1; u ) 

+64H (1, 1, 0; n ) - 32tf (1, 1, 1; n ) + (1, 1, 2; n ) - 12H (1, 2, 0; n ) 

-3H (1, 2, 1; n ) + 48# (2, 0, 0; n ) + 12ff (2, 1, 0; n ) - 9ff (2, 1, 1; n ) 

+48i/ (1, 1, 0, 0; n ) - 24H (1, 1, 1, 0; n ) + 9H (1, 1, 1, 1; n ))) . (5.22) 



6. Check of A° n ~ 

QgQ 

Our new result for the integrated NLO massive quark-antiquark antenna ^Q g Q given in 
eq.(5.10) can be tested with expressions known in the literature. The integrated antenna 
function -^Q g Q can be regarded as the order a s part of the real radiation correction to the 
decay rate of a virtual photon into a massive quark-antiquark pair, 7* — > QQ. By adding 
the real corrections obtained with -^Q g Q to the order a s part of the virtual corrections for 
this decay rate, known in the literature, [39, 40, 41] we are able to reproduce the known 
result for the total hadronic decay width at order a s [41, 42]. The details of the comparison 
are given below. All formulas are restricted to vector coupling only and one heavy flavour 
/ of quarks Q. 

The decay width of a virtual photon (7* ) into a quark-antiquark pair is given by, 



= - a E cm N c Q}r v NS (f). 



with the colour factor Cf given by Cf 



2N C ' 



N c the number of colours and, 



(6.1) 



r V NS (f) = 
K v 

A(v) 



X + Cf—Kv 

if 33 , 22 2 _ . 

- L4(n)+ 24 24 
v \ 1 — v z /3 



(l + v 2 



+3v In 



Li? 



1-v 2 



1 



1 2 



l + v 
vlnv. 




(6.2) 
(6.3) 



(6.4) 



is used as the dimensionless variable and v = J¥q. This result 



In this equation, tq = 1— w- 
is obtained from the decay width of Z boson into a quark-antiquark pair [42] by appropriate 
coupling replacements. It can be reproduced by computing the vertex correction ^* d ver t ex 
and the real radiation correction ^y* d Teal separately. The corresponding formulae will be 
given below. 
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The vertex correction T!l f * d vertex is given in [39, 41] in the form of the QCD vertex 
amplitude V^c 2 corresponding to the decay of a virtual photon with momentum p\ + p 2 
into a quark and an antiquark with colours c\ and c 2 (and a^ u = | [7^, 7"]) 



V£ C2 (pi,P2) = u Cl (p l )T^ iC2 (pi,p 2 )v C2 (p 2 ), 
r^ca(Pl»P2) = -iY / 4^Q / 5 clC2 (V x ^ 



(6.5) 



2m \ m 



(6.6) 



where Q/ is the quark charge in units of charge e = y/Aira for one flavour /. The hadronic 
vertex correction T^ vertex then reads, 



-phad ' ' 

1 7* ,vertex r, -rp 9 



spins, colours 

$2 is the integrated two-particle phase space for two massive particles given by, 

* 2 (m,m) = 2-^(.)-^i^f^(r ) 1 /^(^ m) - 
and the amplitude squared is, 

E KkaW = N c F c 2 m (|F| 2 (2(3 - r - 2e) 

spins, colours 

+ \F 2 \ 2 (-^— (3 - 2r - 2e(l - rO)) + 3? (FF*) (4(3 - 2e)) 
1 - r 

The form factors Fi and F 2 are given up to order a s and e° by, 

F2 = C(e)^ R 



2vr 



with the normalisation factor C(e) given by, 



(6.7) 



(6.8) 



(6.9) 

(6.10) 
(6.11) 



C(e) = (4vr)T(l + e) 



and, 



Fi 



-1+1 



1-y 1 + ; 



#(0;y) 



1 



) 



(6.12) 



2 



#(0;y) 



( X " " ITy) [ 4C(2) " 2i7(0; y) ~ m 0; y) " 2H (1 ' 0; y) ] ) ' (fU3) 
i^(0;y) N 



1-y 1 + 
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We quote the one-loop renormalised vertex corrections for E 2 m > Am 2 given in [39] up to 
order a s and e° respectively. (Imaginary parts are omitted because they do not contribute 
to O (a s )). The dimensionless variable y is defined by 

The real emission contribution to the decay width r^? d real is given by 

r 7Veai = ^\\ E \ M T^\ 2 J ^WW 6 ) (2C F A° QgQ ) , (6.14) 



spins, colours 

while the averaged leading order matrix element squared proportional to lA^, 6 ^^ reads, 
£ \ M r"J 2 = i^Q})N c E 2 cm 2 (3 - r - 2e) . (6.15) 

spins, colours 

By expanding in e our result for Aq^q given in terms of hypergeometric functions we could 
show that, 

-phad _|_ -phad phad /n -i /?\ 

^ 7*, real ~r ^ 7*, vertex 7* ^D.IDJ 

giving us a strong check of our integrated massive final-final QQ antenna A° n ^ . 



7. Conclusions and Outlook 



We have generalized the antenna subtraction method, originally developed for massless final 
states to real radiation off massive final state partons at the next-to-leading order level. The 
main building blocks of the antenna subtraction method are the antenna functions. Those 
encapsulate all unresolved radiation emitted between two colour-ordered hard radiators. 
As such those functions account for all unresolved radiation of the corresponding QCD 
matrix elements. In Section 3, we have presented the final-final antenna functions for the 
following cases: the presence of two massive final state radiators of equal masses or the 
presence of one massive and one massless radiator. Section 4 contained a list of all non- 
vanishing unresolved limits of the antennae functions presented in Section 3. Explicit phase 
space factorisation and parameterization formulae were presented in Section 5 where all 
massive final-final antennae functions were integrated over their corresponding phase space 
measure. 

An important extension of subtraction methods is the combination with parton shower 
algorithms [20, 21, 22, 23], thus allowing for a full partonic event generation to NLO 
accuracy. So far, our massless final-final antenna functions are part of the parton shower 
VINCIA [18]. With the formulation of the antenna subtraction method for final state 
radiation off massive fermions presented here, it will become possible to construct antenna- 
based parton showers involving massive final state particles. 
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